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Abstract 

In this paper we introduce a method for finding a time independent 
Hamiltonian of a given dynamical system by canonoid transformation of 
canonical momenta. We find a condition that the system should satisfy to 
have an equivalent time independent formulation. We study the example 
of damped oscillator and give the new time independent Hamiltonian for 
it, which has the property of tending to the standard Hamiltonian of the 
harmonic oscillator as damping goes to zero. 



1 Intoduction 

It is known that the Hamiltonian description of the system is nonunique. Gen- 
erally, one can change phase-space coordinates (q, p) and in some cases it is 
possible to find a Hamiltonian for these new coordinates. The well known class of 
such transformations are canonical transformations which preserve Hamiltonian 
structure for all possible Hamiltonians (0, But there is also much wider 
class of transformations, called canonoid which preserve at least one Hamilto- 
nian structure. Theory of time canonoid transformations is frequently investi- 
gated in a context of inverse problem of the calculus of variations ([S], 0] and 
references there). General theory of time independent canonoid transformations 
can be found in [3]. 

The main goal of this paper is to develop a method for finding a time indepen- 
dent Hamiltonian for a given system. We restrict ourselves to transformations 
which do not change q, because we could have some physical (e.g. quantum) 
reasons for preserving configurational coordinates. (Without this restriction 
the problem could be trivial, as the transformation to coordinates of integrals 
of motion would yield the Hamiltonian constantly equal to zero). This kind of 
canonoid transformation sometimes is called fouling transformation (0). 
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The paper is organized as follows. First, (section 2) as an introduction, we 
present theory of canonoid transformations preserving q. Then (section 3) we 
find additional relations which the transformation must satisfy to provide the 
transition to phase space of the time independent Hamiltonian. Wc show that 
such a transition is not possible for all Hamiltonians, and we find equation which 
is necessary condition for the system to have a time independent Hamiltonian. 
Finally (section 4), wc test all this methods on the damped oscillator. Wc find 
the time independent Hamiltonian for this system, which tends to the classical 
Hamiltonian of the harmonic oscillator as damping goes to zero. 



2 Canonoid transformations preserving q 

Consider two dynamical systems. One described by coordinates (g, p, t) and the 
Hamiltonian Ho{q,p,t), and the other expressed by coordinates (q, 7r,t) and 
the Hamiltonian H{q,-K,t). Let us also have the mapping p{q,p,t) between 
momenta p and tt, such that {q,p,t) i— » {q,p{q,p,t),t) is diffeomorphism. We 
will call Hq and H equivalent if each solution of Hamilton equations for Hq is 
also (by means of p) solution of Hamilton equations for H and vice versa. One 
can easily check that if both Hamiltonians lead to unique solutions (in terms of 
initial conditions) then it is sufficient only to check the transition from Hq to 
H. Condition of equivalence for the inverse transformation is then automatically 
fuUfiled (at least locally). 

The Hamilton equations for Hq read 



dpi 



q=q(t) 
p=p(t) 



dqi 



q=q{t) 
P=p(t) 



(1) 



Now H is equivalent to Hq if for every solution q{t), p{t) of these equations, 
q{t) and 7r(t) = p{q{t),p{t), t) are solutions of Hamilton equations for H 



dg» , . _ dH_ 

dt ^ ^ ~ dm 



q=q{t) 
7r=7r(t) 



dTT, _ _ dH_ 
dt ^ ' ~ dqi 



q=q(t) 
Tr=7r(t) 



(2) 



Comparing these two sets of equations one gets the following relations 



OHq 



dpi 



q=q(t) 
P=P{t) 



dH 



q=q{t) 

n=piqit),p{t),t) 



and 



{Pj,-f^o}|q=q(t) + ^ 



dpj 

dt 



q=q(t) 
p=p(t) 



p=p{t) 

Then, in terms of the coordinates, we obtain 



dH 
dqj 



q=q(t) 

■^=p{q(.t),p{t),t) 



dHQ 



dpi 



q=q 
p=p 



dH 
d-Ki 



q=q 

7r=p(q,p,t) 



(3) 



2 



and 



{P],HQ}\q=q - ^ 

p=p at 



q=q 



dH 



q=q 

■^ = P{Q:P,t} 



(4) 



These are our fundamental relations and we can treat them as the equations for 
the new Hamiltonian H and for the mapping p. We can easily transform them 
to the following form 



dHo 



dpi 



dpi 



d 



<i=q dpj dp'j 
p—p 



H\q = q 

■^=P{<l:P-t) 



and similarly 



dpj 

{Pj,Ho}\q=q~ 



q=q 
P=P 



dpi 



dp^ d 

Tr~ = Tr~ I H\q=q 
q=q oq. oq. 
p—p 



■^=P{Q,P,t) 



(5) 



(6) 



Above equations can be viewed as definitions for the derivatives of H\q=q 

■^=p{q,p,t) 

But such definition could make sense if the integrability conditions 



d^ ( , 
dqjOqk V n=p(q.p,t] 



92 



92 



dpjdpk 



52 



dpjdqk 



H\q=q 

T^ = Pi<l,P:t) 

H\q = q 

■^=p{q,p,t] 



dqudqj 

92 
dpkdpj 

d^ 
dqkdpj 



H\q=q 

■^ = P{l:P:t) 

H\q = q 

■^ = Piq:P:t) 

H\q=q 

n=p{q,p.i) 



are fulfilled. This yields the following equations 



_d_ 

dqk 



d 
dpk 

dH 

dpr 



{Pi,Ho} 



{pj,Ho} 



dPj_ 
dt 



dt 



_d_ 

dqj 



{Pk,Ho} 



= 0, 
= 0, 
= 



dpk \ 
dt ) 



where qm{q,P, t) — and we use the Lagrange bracket defined by 



d(j)m dlpm d(j)r, 



(7a) 
(7b) 
(7c) 

(8a) 
(8b) 
(8c) 

(9) 



with summation over m. Let us remark that for one degree of freedom (one q 
and one p) equations (|8a|l - l|5c|) reduce to one relation 



dp 
op 



d^p 
dtdp 



0. 



(10) 



Equations (|8a|l - l|5c|) do not contain H and are linear with respect to p, so they 
are remarkably simpler then and l@J. 
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Once (|8a|l - (jSc)l are solved with respect to p we can use the solution to find 
H\q=q . Using Poincare lemma to the form defined by right sides of ((Sj 

and © we can write 

= [ (-{pj,Ha} - ^ + ^It^) |q=qo+r(q-qo) dr {q^ - qoj) 

P=PO+T{p-po) 



7v=piq,p,t) Jo \ dt dp, dq 

'dHo dp. 



dpi dpj 



q=qo+r(q-qo) [pj ~ POj) + F{t), 
P=Po+t(p-po) 

(11) 



where (qo,Po) is some arbitrary point for which assumptions of Poincare lemma 
hold (the form we integrate must be well defined in a star-shaped neighborhood 
of (goiPo) for every moment of time t). 

To summarize this section let us write down a general scheme for finding 
the new Hamiltonian equivalent to Hq, expressed by the same configurational 
coordinates q and some new momenta tt. 

• First one must solve the equations H8a |l - (|Sc|l to find p i.e. possible trans- 
formation of momenta. 

• Using formula one should find H\q=q 

■^=p{q,p,t) 

• Finally, one must find the inverse function p~^ and express the new Hamil- 
tonian H as a function of (q, tt, t). 

The equations presented in this section are analogous to these in where 
canonoid transformations with dependence in q (but not in t) are studied. 

3 Time independent Hamiltonian 

We are going to find additional conditions under which the new Hamiltonian H 
appears to be time independent. First note that from equality 

Pj~^{Q,PiQ,P,t),t) ^ Pj. (12) 

one can obtain two useful relations. The first one is the well known formula for 
the derivatives of an inverse function 



dpT' 



q=q 

■^=p{q,p-t) 



where (f) "Ss the inverse matrix of the invertible matrix with elements a., = 
The second relation is of the form 



dp 



dt 



1=1, \dpjj^ dt 



-^^piq-.p-t) 



dp\ ^ dpn 



(14) 
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and can be derived by differentiation of (|12|l with respect to t. Consider now 
the equation If we compound it with we obtain to be 



dH 



q=q 



ri2 TT 

But we want our new Hamiltonian to be time independent, so — 0, and 

consequently 



q=q 



dt 



dpidt 

Then we compound above expression with p 



q=q 

p=p-^{q,-K,t} 



0. 



dt 



q=q 

■K=p{q,p,t) 



dpidpj 
After using 114(1 we get 

d^Ho ( dp\~^ dp. 



dpidpj \dp 



dt 



dpidt 



dpidt 



Let K be the inverse matrix of the matrix 



exists unless 



there are constraints imposed on velocities q). Then we can rewrite derived 
formula in the final form 



dpi _ dpi_^ d^Hp 
dt dps dpidt 



(15) 



This is a linear equation for p that must be fulfilled together with (|8a|l - jSc|) to 
lead to a time independent Hamiltonian. 

Analyze equation Q and perform same procedure as before. First we come 
over to 77, then we differentiate the result with respect to t and finally we come 
back to original momenta p. The result reads 



dP] 

dt 



Ho 



Pj> 



Pj 



dpr, 



dHo\ 


' dt 


7" 


dpVn 




dt 


q=q 



dPj 

dpm 



-K=p(q.p.,t) 



Ho 



dp. 



dt 



q=q 

d^Pj 9p-^ 



+ 



d% 

dt'^ dpmdt dt 



q=q 

■^=p{<l,P,t) 



Substituting into (d)) we get 
dp^ 



dt 



q=q 

n=p[q,p,t) 



d^Ho 
' dpvdt 
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This formula and (|14|l lead us to the following expression 

Opm OpvOt 



1 



Jo 




p.^ I " _ ( "^J K„iv " I + " ^'^ Kmv " = 0. 

^'dpjnj "^^ dpi,dt dt \dpm dpi,dt J dp„idt "^^ dpvdt 



2 3 3 

Now, after cancellations in groups indicated by numbered braces, we may write 
dp, ( d^Hp \ dp~j cFHo dp, JPHo^ d^Ho , 



OPm [ OpyOt ) Opm OqmOt Opm Oq„iOPs OpyOt 

' , ' " « ' 

1 2 



dpm dt \ dpydt 

" V ' 

3 

Finally, elements of invertible matrix, we get the main result 



dqidt dqidps dpvdt [ dpvdt J dt \ dpvdt 

In this equation there is no H, p or tt. Consequently we interpret it as a 
condition that Hq should satisfy to have a time independent counterpart. 



4 Examples 

4.1 Pi 



^ dqi 



— is a solution to H8a|l - for arbitrary 



One can easily check that Pi ~ Pi -r -g^ 
f{q,t) and every Hq. Integrating (fTT|l around {q,p,t) we get 



df df 
H\q=<, =Ho{q,p,t)~^-Ho{0,0,t)+ ^ 

7T—p(q,p,t) OZ Oi 



Fit). 



q=0 



Then expressing new Hamiltonian as a function of tt and putting F{t) 
Ho{0, 0, i) - If we obtain 



q=0 

H{q,7v,t) = i/o(g,7r- - ^. 

One quickly shows that such transformation of Hamiltonian corresponds to the 
standard ambiguity in Lagrangian for a system (namely adding term §^qi + ^ 
to Lagrangian). Notice that this transformation is canonical because it holds 
for every Hamiltonian. 
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4.2 The damped oscillator 

The damped osciUator is a commonly used system for testing theoretical ideas 
(see e.g. [7j - covering topics closely related to ours, or IB] and references there). 
Let us find a time independent Hamiltonian for it. We additionally want our 
new Hamiltonian to become the standard Hamiltonian of the harmonic oscillator 

( 2 2 2 \ 

^ + ^^^-^ 1 if there is no damping. There exists an example, given by Havas 

of the time independent Hamiltonian for damped oscillator. But the Havas 
Hamiltonian yields quite strange form in the case when j3 — 0. As a starting 
point we use the Hamiltonian 

^. = ^[e-'^V^ + (/3'+c.W], (17) 

which is time dependent but for /3 = is exactly the Hamiltonian of the har- 
monic oscillator. After simple calculation we observe that ()17|) satisfies condition 
P6() . In our case the equation (|15|l becomes 

^ OR ^P n 
— +2(}ps— = 0. 

at ops 
General solution to this equation is 

p^v (g, z) , (18) 

where z = PsS"^'^* and v{q,z) is an arbitrary (sufficiently smooth) function of 
two variables. But p must also be a solution to ifTUI) . If we put (fTH|) to ((117)) we 
obtain after some transformations the following relation 

^^z-^^{{P'+u;')q + 2Pz)=2Pg, (19) 
where g{q, z) — ^ . The general solution to this equation is given by 



1 



q'^uj^ + {z + j3qY ' 



x^(^_i===exp('^arctanfi±^n I, (20) 



\Jq^^'^ + (2: + (iqf \ ^q 

where ip is an arbitrary function. To prove this one observes that ti — g{q^bj'^ + 
(z + l3qY ) and t2 — , ^ exp ( - arc tan ( ) ) are two independent 

integrals of the system of equations 

dq dz dg 

~ ~ ~ (/32 + u;2)g + 2l3z ~ 2^' 
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From ij^OI) and (dJ we have 

dp g-^l3t 



2 



^ , , , . , . . exp - arctan . (21) 



But for /3 = the Hamiltonian should remain unchanged so without damping 
one should have p{q,ps, t) — ps and — 1. Let us write H21|l for P — 0: 

dp 1 f 1 \ 

-■0 



dps q^uP- +p1 y ^g2^2 p2 J ■ 

This allows us to choose = It is not the only possible choice as ip can 
also be a function of (3 and above equation gives us information only for P = 0. 
Using this form of ip we obtain 

dp _2;3t _____ f 2/3 fp,e-^0* + Pq 



e exp arctan — — . (22) 

OPs \ \ ujq J J 

In what follows we use the following notation: 

• B{a, b) ^ t°-i(l - t)''-^dt is the Euler beta function, 

• B{z, a, b) = <°^^(1 — tf'^^dt is the incomplete beta function, 

• I{z, a, b) — B{z, a, b)/ B{a, b) is the regularized incomplete beta function, 

• I~'^{z,a,h) is the inverse of the regularized incomplete beta function, i.e. 
function providing solution of equation z = I{s, a, 6), with respect to s, 
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Using the relation arctan(x) = ^ilog ^f^f^ and integrating we get 
Piq.Ps,t) = 2iujqB{T-,T+)x 

In this result we chose integration constant taking p{q, 0, t) = 0. Note that p 
is real. The inverse transformation can be written as follows (we abandon tt 
as a symbol for new momentum because it could be confusing in context of 
trigonometric functions; instead we use p„). 



P 



-l(q,p„,t) = e2^*X 



X 2iujq 



/-I ( ^ +l(lr-,T-,T+] ,r-,r+ ) - - 

' 2iujqB{T-,T+) V2 y / 2 



-I3q 
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Applying the formula one finds 

H\q=q = ie 

n=p{q,p,,t) 2 dp. 



The form we have integrated is not well defined in (0, 0, t) but one can check that 
the above result holds true. ^ 
Let us pass to momentum p„ 



the above result holds true. We left term to shorten further calculations. 



2 



-2/3t 

dps 
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q uj X 



X 1 



2r 



'( 2^c.g4r-,r+) +^r 



But ^ 
(see M) 



dI-Hz,a,b) 
dz 



1/ 



dp 



By using the formula for derivative of / ^ 



i3(a,6)(l-/-i(z,a,5))i-''(/-i(0,a,6)) 



l-a 



(23) 



we obtain 
dp 

dpn 

Finally 



2iujqB{T ,T+) 



- 1 



H{q,Pn) = 



1 - r 



, T f ^ - - + \ - + 

2iujqB{T-,T+) V2 



2iu!qB{T , 



(24) 



One can check by direct calculations that our new Hamiltonian is a real function 
of q and Pn and it gives equation of motion of the damped oscillator. For P — 
(and consequently t+ = = 1) (|24|l reduces to ^ (pf^ + (7^cj^). Notice that 
there is a discontinuity in 1241) for q = 0. 



5 Conclusions 

We have developed method for generating new Hamiltonians by changing canon- 
ical momenta. We have applied this method to find the condition H16() for exis- 
tence of a time independent Hamiltonian and to show how to construct it. The 
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ambiguity in Hamiltonians is especially intriguing in quantum context. We have 
given example of the damped oscillator and have found the new Hamiltonian 
for it. We believe that this new Hamiltonian can yield some interesting results 
at quantum level. 
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